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Abstract: We calculate the baryon asymmetry of the Universe which 
would arise during a first order electroweak phase transition due to minimal 
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baryonic excess, for resonable KM parameters and mt in the expected range, 
and plausible values of bubble velocity and other high temperature effects. 
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The existance in the present-day universe of an excess of matter over 
antimatter has long been recognized to be a fundamental problem in cos- 
mology. A great deal of theoretical work has been done, and a number of 
possible explanations have been advanced. So far, most of these proposals 
have required physics beyond the standard model. Indeed, the existance of 
the baryon asymmetry of the universe (BAU) has widely been considered 
one of the most compelling pieces of evidence that the standard model is 
incomplete. 

In this Letter we consider the possible production of a baryon number 
excess as a result of purely minimal standard model (MSM) processes, dur- 
ing the electroweak phase transition. Our calculation is realistic enough to 
make it clear that the baryon asymmetry arising from this mechanism can 
be responsible for the observed baryon density to entropy ratio, hb/s ~ 
(4 — 10) lO^"*^^, for values of the top mass and KM parameters in the cur- 
rently favored ranges. In order to settle definitively whether this is the actual 
origin of the baryonic asymmetry existing today, it will be necessary to im- 
prove our quantitative understanding of the ew phase transition and of the 
behavior of quasiparticles in the high temperature plasma, as well as refine 
and extend the analysis reported here. 

For an excess of baryons to develop in a Universe which initially has zero 
baryon number, the following conditions, first enunciated by Sakharov, must 
be met: 

1. Some interaction of elementary particles must violate baryon number, 
since the net baryon number of the universe must change over time. 

2. C and CP must be violated in order that there is not a perfect equality 
between rates of AB ^ processes, since otherwise no asymmetry could 
evolve from an initially symmetric state. 

3. A departure from thermal equillibrium must play an essential role, since 
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otherwise CPT would assure compensation between processes increas- 
ing and decreasing the baryon number. 

We briefly summarize several features of the standard model which are nec- 
essary to understanding how these requirements will be met: 

1. In the standard model, quarks get their masses as a result of their 
Yukawa couplings to the Higgs field. When the Higgs field has a non- 
zero vacuum expectation value (vev), quark masses are induced which 
are proportional to their couplings to the Higgs field, times its vev. 
There are off-diagonal Yukawa interactions, in which quarks of dif- 
ferent generations couple to one another through the Higgs field. In 
general the couplings are complex, and for three generations there is 
one physically-significant phase, 5cp- A non-zero value of 5cp is imag- 
ined to be responsible for the CP violation observed in the kaon system. 
The KM matrix describes the mixing between generations and contains 
the phase of the Yukawa couplings. 

2. With three generations the phase in the KM matrix could be rotated 
away if any pair of quarks of the same charge were degenerate in mass, 
or if any of the mixing angles vanished. Thus CP violating effects are 
significant in particle physics only when a relevant scale is small enough 
to be of the order of magnitude of the splitting in mass between, e.g., 
iris and m^. The K° system is an example of this, where the lack of 
degeneracy between d and s is evidenced by mx 7^ w^tt- Moreover CP 
violation in the MSM vanishes together with 

J = sin{6i2)sin{6i:i)sin{623)sin{6cp), (1) 
using the particle data group convention for the KM mixing angles. 

3. Although the standard model Lagrangian conserves baryon number, 
quantum effects produce an anomaly which leads to baryon number 
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violating transitions. While the rate of these transitions is negligible 
at T = 0, at high temperture their rate ~ exp[— '^'^^') ] [g| . Thus 
above the electroweak phase transition the baryon-number-violating 
transition rate is rapid compared to the expansion rate of the universe. 
We will require that the vev in the low temperature phase is large 
enough that baryon violation is "turned off", allowing the asymmetry 
which has been produced during the transition to survive. In the MSM 
where the only undetermined parameter of the Higgs sector is the Higgs 
mass, this requirement may lead (if non-perturbative thermal effects are 
unimportant) to an upper limit on the Higgs mass.0 

At temperatures above the electroweak phase transition, the vev of the 
Higgs field vanishes it takes on a constant, non- vanishing value in the low 
temperature phase. We require the phase transition to be first order, al- 
though this depends on the Higgs mass and is not certain to be true. During 
the phase transition, "bubbles" of the new < $ >7^ phase nucleate and 
expand to fill the universe. This departure from thermal equillibrium satis- 
fies the third Sakharov requirement. Bubble properties, such as their wall 
thickness or their expansion velocity, are not well understood at present. A 
number of phenomena could produce a baryon excess during the ew phase 
transition. For definiteness we concentrate on a specific mechanism; see 
for others. 

As a result of their thermal motion, quarks and antiquarks in the neigh- 
borhood of the bubble wall propagate through it. Since their masses result 
from their interaction with the vev, they see the bubble wall as a potential 
barrier and scatter from it. We model this process in detail^, keeping the 
plasma masses of the quarks and antiquarks which originate from their inter- 
actions with the gauge and Higgs particles present in the heat bath, as well 
as treating quantum mechanically the process of their interaction with the 
bubble wall of Higgs field. Due to the spatial variation of the effective CP 
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violating phase, which comes about because the physical eigenstates depend 
on the interplay between flavor dependent thermal effects and the masses 
induced by the changing Higgs vev, there can be a difference between the 
reflection and transmission coefficients of quarks and antiquarks. We report 
below on our computation of this asymmetry, in the one-dimensional prob- 
lem which results when quark motion parallel to the bubble wall is ignored. 
Due to the fact that Lorentz invariance is broken by the thermal medium, 
momentum components parallel to the wall could be dynamically important, 
but that issue will be left for future work. 

The total baryonic current is conserved in quark scattering with the bub- 
ble wall. Nonetheless, if there is an asymmetry in reflection and transmis- 
sion coefficients the wall would separate particles and antiparticles, with, 
e.g., quarks flowing preferentially toward the low-temperature phase and an- 
tiparticles toward the high-temperature side. In the high temperature phase, 
sphaleron transitions operate to equilibrate the antiquark excess 0, convert- 
ing most of them to quarks and leptons. But as long as the vev in the low 
temperature phase is sufficiently large, the sphaleron transition rate is too 
low to keep up with the expansion rate of the universe, and the quark excess 
is preserved until now.|^ 

As noted above, if two same-sign-quarks are degenerate in mass there 
is no CP violation, since the phase 5c p can be removed from the KM ma- 
trix by a change in the definition of the overall phases of quark fields. This 
fact manifests itself in the present context by a tendency for different flavor 
eigenstates to have canceling contributions to the baryonic asymmetry cur- 

^ The idea that the BAU could resuh from the separation of a quantum number by the 
bubble wall combined with equilibration of this quantum number in the high temperature 
phase due to baryon-number- violating sphaleron processes, originated with Cohen, Kaplan 
and Nelson, employing a lepton-number violating interaction with the bubble wall|^. The 
idea that MSM interactions of quarks and antiquarks with the bubble wall could directly 
cause a separation of baryon number is due to MS0], where the element of including 
thermal effects is also introduced. 
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rent. For instance, if the reflection probability at a given energy for a di is 
greater than that for an anti-di, the reflection probability of an sl will be less 
than that for an sl by a nearly identical amount. In the limit — — 
the compensation is perfect, when reflection involving 6's is included. An 
estimate of the residual CP violating asymmetry for typical quark energy 
~ T is (m^ — m^)(mj — m^)(m^ — ml) {ml — m1){ml — ml){m1 — m^/T^"^ 
times J, ~ 10^^^ 0. For this reason it has commonly been believed that the 
MSM cannot be responsible for the BAU. 

The new observation of the present work is that the important quark 
energies are not ~ T, but rather energies such that the s quark is totally 
reflected but the d quark is partly reflected and partly transmitted, max- 
imizing the dynamical difference in their contributions to the asymmetry 
current. Taking into account thermal effects on the quark propagation (see 
below and ref. 0), one finds that the pt for which s quarks are reflected 
is Pt ~ ^^^Jy^T ~ 10~^T. The range of pt with s but not d reflection is 



The net baryonic flux through the wall is proportional to the group veloc- 
ity of the quasi-particles perpendicular to the wall (~ |) and to the asymme- 
try in Fermi-distributions on the two sides of the wall due to its motion with 
respect to the themal medium (~ 2v^^ ~ |f -y/^^ for small wall velocity,^, 
where u = uj{p) is the quasiparticle dispersion relation). The strength of CP 
violation in this region is just J, given in eqn ([^). Thus putting all the factors 
together we expect 

2 l2TTas f ms — md\ Jf , . 

where Nf,ff ~ 100 is the total number of degrees of freedom. Global fits 
to determine KM parameters place J in the rangef^ (1.4 — 5.0)10"^, so we 
expect rih/ s ~ (2 — 8)10^^^f , which can be of the right order-of-magnitude 
to account for the BAU. 
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Having outlined the way in which CP violation in the quark mixing ma- 
trix can lead to a present-day baryon asymmetry, we next calculate it more 
quantitatively, in the quark-antiquark separation mechanism. There are two 
essential effects to be included: The interactions of the quarks and antiquarks 
with the plasma of gluons, ew gauge bosons and Higges, and the quantum 
mechanical scattering of the quarks and antiquarks from the bubble wall in 
the Higgs vev. The effects of the interactions of the quarks and antiquarks 
with the gauge and Higgs particles in the heat bath are most efficiently taken 
into account by changing variables to a quasi-particle description. The prop- 
agators of the quasi-particles have been determined to one-loop accuracy, 



neglecting internal masses, by KlimovQ and Weldon|T^. The quasiparticles 
have interesting and unfamiliar behavior, but space limitations prevent us 
from describing them further here. These and other details are given in a 
longer article [^]. Here we record only the equation of motion which these 
quasiparticles obey in the wall rest frame, in the limit which is relevant to 
total reflection of the strange quark, momentum <^ tu, treating only the mo- 
tion normal to the bubble wall, and in the limit of small plasma velocity with 
respect to the wall: 

u{l + aL + I3l) + + ai) KMa{x) \{L\=n 

(3) 

where aL,R, Pl,r and Md{x) are 3x3 matrices and L and R are 3 component 
spinors in flavor-space. M^ix) is the Higgs-induced mass at T^. For charge 
-1/3 quarks in the unbroken phase, in the gauge basis where the interactions 
of quarks with the W and Z are diagonal: 

An as 47r Km\K'^ sin^ 6^^ 

" ^3 6 3 32 ^Ml 9 ""'^u^' ^ ^ 

where and are diagonal matrices of the charge -1-2/3 and -1/3 masses 
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at T = 0, and K is the KM matrix. For momenta small compared to 
Pl,r = '^OiL,R- More general expressions can be found in ||^. 

We can solve these equations analytically in two limits: no mixing and 
zero wall thickness with small mixing. For more realistic cases they must 
be solved numerically, although having the exact cases to verify the correct- 
ness and accuracy of our numerical solutions is very useful. Details of the 
analytical and numerical results are given in ref. |^. We find that when 
the energy is such that neither or both d and s quarks are totally reflected, 
the difference between reflection probabilities of quarks and antiquarks, after 
summing over all three flavors, is extremely small. (Less than our numerical 
integration accuracy of one part in 10^°.) However, as expected on the basis 
of the heuristic discussion above, the asymmetry is substantial in the narrow 
energy window in which the s quark is totally reflected but the d is not. 

The figure shows A, the difference in the reflection probabilities for right 
chiral quarks and antiquarks incident from the unbroken phase, summed 
over flavors, in the interesting range of energies. We take as representative 
values rrit = ISOGeV, nic = 1.6GeV, m„ = .OOSGeV, nib = 5GeV, = 
.15GeV, rrid = .OlGeV, S12 = .22, S23 = .05, S13 = .007 and sin{Scp) = 1. 
For the calculation of this figure the wall velocity, v, was zero and the wall 
thickness was 10/Tc, a popular value. Taking the wall to be narrower does not 
significantly change the result; taking it a factor of three thicker increases the 
result by a factor of two. This dependence on wall thickness is not surprising: 
even in the thin wall limit there is a non-trivial CP conserving phase shift to 
interfere with the KM CP violating phase. The asymmetry in the reflection 
probabilities increases when the effect of the flow of the thermal medium is 
included|P|, so that A is almost a factor of five larger for v = .25 than it is 
for the V = case shown in the figure. 

The upper pair of peaks occupy the energy range in which the strange 
quark is totally reflected. Note that the width in energy of this region is ~ .1 
GeV, just the mass of the strange quark at that temperature. The "notch" in 
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the middle, of width ~ .006 ~ md{T), is the region in which the down quark 
is also totally reflected and GIM cancellation is perfect, as expected. The 
unfamiliar feature that total reflection occurs for a range of energies, rather 
than for all energies less than some value, results from the unusual properties 
of the quasi-particle dispersion relation^, but is not essential to our result. 
Our analytic calculation in thin- wall, small mixing approximation]^ provides 
an adequate description of reflection in this region. At lower energy there 
is another region of a different character, involving level crossing between d 
or s and 6r. It would not be present if mixing were absent. Its width is 
~ mb(T) sin(6'23). 

We have checked that the asymmetry vanishes as pairs of masses are 
brought together. When — >• it arises by the squeezing away of the 
width of the upper peaks and the diminution of the magnitude of the lower 
peak. When masses in the charge +2/3 sector are brought together, or nib 
rrid^s, the magnitudes of the peaks decrease appropriately. We checked that 
the result vanishes as mixing angles are taken to zero, although in the physical 
range of 623 the result is non-linear, increasing by 40% as 623 is changed 
from 0.05 to 0.06, and changing sign for ^23^0.03. The dependence is 
interesting: for low values, mt<110GeV, the integrated asymmetry has the 
opposite sign as for m^^llOGeV. It reaches its maximum value for nit ~ 
210GeV, where it is more than four times greater than for mt = ISOGeV, 
then decreases for larger mj. For a more detailed discussion and additional 
figures see ref. P]. 

The net L baryonic current resulting from the asymmetry in reflection 
coefficients for i?'s incident from the unbroken phase and L's incident from 
the broken phase is obtained from A as follows 0: 

Given the net baryonic current flowing through the bubble wall, correspond- 
ing to a preferential flow of baryons in one direction and antibaryons in 
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the other, we next wish to determine the resuhant ub/ assuming sphaleron 
transitions operate on L chiral quarks and antiquarks to balance the chemical 
potentials in the high temperature phase, but are completely suppressed in 
the low temperature phase. To make an accurate estimate requires a detailed 
study of the non-equillibrium statistical physics of the problem, however we 
can obtain a conservative estimate of n^/s as follows. Suppose that the wall 
velocity is low, so that diffusion permits a back-current of baryon number 
to be established, which acts to replace the anti-baryon number which is 
being destroyed in the high temperature phase at a rate F by the sphaleron 
transitions. For sufficiently low velocity the problem is essentially static and 
the result is = 3J^p/(p), where p = /(p) = 1 for p ^ 1 and 

/(p) = |p for p ^ 1. The physical importance of p is clear since the typical 
distance from the wall, of a particle scattered at t = into the unbroken 
phase, is Voi — vt. Thus ^ is the typical time in which the sphaleron 
transitions can act on that particle before it is enveloped by the expanding 
low-temperature phase. The correct values to take for the diffusion length, 
D, and sphaleron rate and wall velocities are very uncertain, but can plau- 
sibly be such that /(p) ~ 1, though a suppression as large as 10~^ is also 
possible [0. 

Since we have computed the current in one-dimension, we divide by the 
one-dimensional entropy for the MSM particle content at the ew phase tran- 
sition: = -^^j^- Taking v = .25 values for reflection probabilities and 
the Boltzman factors, we find 

hb/s ~ 4 10-'^ ( /(p), (5) 

' V.22 X .05 X .007/ ^ ' 

for TJit = ISOGeV, 623 = 0.05, and inverse wall thickness lOGeV. (Until 
the result of the full three-dimensional calculation is known, one should take 
this result as an order-of-magnitude estimate.) Since the prediction increases 
rapidly for larger values of rrit and wall velocity, there seems to be ample mar- 
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gin within the favored ranges of these quantities^ to tolerate some suppres- 
sion from the uncertain overall factor sin 5c p ( '^'007'' ) / (p) ^^nd the difference 
between the one- and three-dimensional cases.Q 

While the sign of sin 6cp is not at present unambigously determined fTl 



a positive sign is favored. In this case (H) correctly predicts a baryon, not 
antibaryon, excess. Within the framework of this model, effects which we 
have not included seem capable only of affecting the magnitude and not the 
sign of this prediction (see ref. 0). Changes in v and wall thickness, and 
changes of the poorly known rrit, 623, and ^13 within their favored ranges, do 
not change the sign of (^). Thus refinements in the treatment of this problem 
may not modify the conclusion, that minimal standard model interactions 
can be responsible for the BAU. 

The most crucial outstanding problems are those associated with our 
still-primitive ability to deal with the high temperature environment during 
the electroweak phase transition: sphaleron rate, wall velocity, quasiparticle 
scattering length are obvious examples. In addition, our calculation of the 
quasi-particle reflection coefficients requires knowing their equation of mo- 
tion in the plasma, which we determined by using propagators calculated 
to 1-loop accuracy in high temperature perturbation theory. Until plasma 
screening and "confinement" effects are better understood, these propagators 
must considered as approximate. Also, a more complete calculation would 
include incoherent as well as coherent processes occuring during the quantum 
mechanical scattering from the vev. 

To summarize, we have argued that already-known physics of the minimal 
standard model may explain the observed baryon asymmetry of the universe. 
A quantitative calculation in a specific mechanism gives the correct sign and 
'^Recent estimates place .1 < v < 1, see § for refs. 

^Of course if there were a fourth generation with a comparable KM structure then the 
degeneracy between d and s would be irrelevant and would be replaced by the degeneracy 
between b and s. Using the analytic thin wall expression!^ yields an enhancement factor 

_ /rrv)4/m^)2/m^)2 3,3324^34 _ iq3 f _ _ SOOGeV. 
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magnitude. The essential new ingredient is not overlooking those regimes of 
quark momenta in which the most degenerate pair of quarks have very differ- 
ent dynamical behavior. If this is the explanation for the baryon asymmetry 
of the universe, then future precision comparisons between observation and 
theory will provide a powerful test of our understanding of the ew phase tran- 
sition, as well as constrain the KM matrix and the masses and generation 
content of the standard model. 
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Figure Caption: A, the asymmetry in the reflection probabilities of 
right-chiral quarks and antiquarks incident from the unbroken phase, for 
zero wall velocity. 
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